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Abstract 

If spacetime has more than four dimensions, ultra- high energy cosmic rays may create micro- 
scopic black holes. Black holes created by cosmic neutrinos in the Earth will evaporate, and the 
resulting hadronic showers, muons, and taus may be detected in neutrino telescopes below the 
Earth's surface. We simulate such events in detail and consider black hole cross sections with and 
without an exponential suppression factor. We find observable rates in both cases: for conservative 
cosmogenic neutrino fluxes, several black hole events per year are observable at the IceCube detec- 
tor; for fluxes at the Waxman-Bahcall bound, tens of events per year are possible. We also present 
zenith angle and energy distributions for all three channels. The ability of neutrino telescopes 
to differentiate hadrons, muons, and possibly taus, and to measure these distributions provides 
a unique opportunity to identify black holes, to experimentally constrain the form of black hole 
production cross sections, and to study Hawking evaporation. 
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I. INTRODUCTION 



The possibility that we live m D = A + n > A spacetime dimensions has profound 
implications. In particular, if gravity propagates in these extra dimensions, the fundamental 
Planck scale Mr> at which gravity becomes comparable in strength to other forces may be far 
below M 4 ~ 10 19 GeV, leading to a host of potential signatures for high energy physics |[[. 

Among the most striking consequences of low-scale gravity is the possibility of black 
hole creation in high-energy particle collisions 0, |, ||, §, || [7|, [jj. In most gravitational 
processes, such as those involving graviton emission and exchange, analyses rely on a per- 
turbative description that breaks down for center-of-mass energies of Mp and above. In 
contrast, black hole properties are best understood for energies above Mp, where semi- 
classical and thermodynamic descriptions become increasingly valid [TIJ. In principle, then, 
black holes provide a robust probe of extra dimensions and low-scale gravity, as long as 
particle collisions with center-of-mass energies above Mp ~ 1 TeV are available. 

Nature provides interactions at the necessary energies in the form of cosmic rays with 
energies above 10 10 GeV. In collisions with nucleons, these cosmic rays probe center-of-mass 
energies exceeding 100 TeV, beyond both current man-made colliders and those of the fore- 
seeable future. Cosmic neutrinos may create black holes deep in the Earth's atmosphere, 
resulting in spectacular signals of giant air showers in ground arrays and fluorescence detec- 
tors 0, [L2[ [13], [14|, jlSf . With a handful of events, standard model (SM) and most alternative 



explanations may be excluded [11], |15| by comparison with rates for Earth-skimming neutri- 
nos [16|, [L7|, [18| , and with more events, black holes may be identified through their shower 
characteristics fL2| . Bounds have been derived [|14], pl| from the absence of such showers in 
current data from Fly's Eye |19j and AG ASA |2(J. For conservative fluxes and the geometric 
black hole cross section, the AGASA data require Mb > 1.3 — 1.8 TeV for n > 4, the most 
stringent constraint to date [15|. At the same time, the Auger Observatory, scheduled for 

may observe tens of black hole events per year JTT], [14], [I3| . Related 

" El- 



completion by 2004 

phenomena related to p-brane production may also be observed in cosmic rays |2T 

Here we examine the possibility of detecting and studying black holes produced by cosmic 
neutrinos in neutrino telescopes. Several large-scale neutrino telescope projects are under- 
way, including IceCube pj in the Antarctic ice, and ANTARES [E3] and NESTOR E5 



in the Antarctic ice, and ANTARES @] and NESTOR 
in the Mediterranean. Among the many possible black hole signatures, such detectors are 
most sensitive to contained hadronic showers and through-going muons and taus from the 
evaporation of black holes produced in the Earth's crust. The possibility of black hole de- 
tection in neutrino telescopes has recently been studied in Ref. p6| ; where possible, we will 



compare our results to those of Ref. ^6f below. For a preliminary study, see also Ref. [^7]] . 

For TeV-scale gravity and conservative flux assumptions, we find that IceCube could 
detect several black holes per year. These rates may be enhanced by larger fluxes, and 
observable rates are possible even given a postulated exponential suppression factor in the 
black hole cross section [^, [2~9fl . The relative event rates in the three channels may differ 
from the SM, and the energy and angle distributions of black hole events are also distinctive. 
These will not only help identify black holes, but may also constrain parameters such as n 
and Mp, and determine if suppression factors in the cross section are present or absent. The 
search for black holes in neutrino telescopes therefore complements black hole searches in 
other cosmic ray detectors, as well as searches for the effects of perturbative gravity processes 
at center-of-mass energies below Mo |50[ [H], [32], |33| . 
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II. BLACK HOLE PRODUCTION AND EVAPORATION 



In D = 4 + n dimensions, gravity is described by the Einstein action 



8ttG 
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-gin 



(1) 



where Go is the D- dimensional Newton's constant. To define the fundamental Planck scale 
Mp, we adopt the convention 

i M 2+n 

1 = . (2) 

8nG D (2tt)« 1 ; 

In the most straightforward scenarios [Q with flat extra dimensions of equal length, TeV- 
scale gravity is excluded for n = 1 by solar system tests of Newtonian gravity. Astrophysical 
bounds |34j] on supernova cooling and neutron star heating provide the most stringent bounds 



for n = 2 (M D > 600 TeV) and n 
bounds are from collider searches for perturbative graviton effects |35 
ray bounds on black hole production 



3 {Md > 10 TeV). For n > 4, the most stringent 

and cosmic 

fl4] , [15| ) |19| ) p0|| . These constraints are each subject to 
a variety of theoretical assumptions (for a comparison and discussion, see, e.g., Ref. [|15|]), 
but the most stringent of these require roughly M D > 1 TeV. 

To determine event rates for neutrino telescopes, we must first model black hole produc- 
tion and evaporation. For production by cosmic neutrinos, we follow the analysis of Ref. ||TT] . 



Black holes produced in parton collisions are typically far smaller than the length scales of 
the extra dimensions. These black holes are then well-approximated by (4+n)-dimensional 
solutions. The Schwarzschild radius for a (4+n)-dimensional black hole with mass Mbh and 
vanishing charge and angular momentum is 
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We assume that two partons i and j with center-of-mass energy \/§ form a black hole of 
mass Mbh = when they pass within a distance r s (s), leading to a geometric cross section 

of si a 

a( i3 ^BE)(s)=7rr 2 s (s) . (4) 



Evidence from analyses of axisymmetric |39[ and off-axis classical collisions, an analysis 
in a simple model framework Q , and a string calculation [42] suggests that this picture 
is valid semi-classically and is not subject to large corrections |43|. Modifications for non- 



vanishing angular momentum and spinning black holes have also been found to be small [15 



However, Voloshin has argued p8, that the cross section could be suppressed by the factor 



where the action is 



S 



47tM B h r s 



(5) 



n + 1 (n+l)(n + 2) ' 

with S the black hole entropy. This implies vanishing cross sections in the classical limit, 
contrary to the evidence noted above. To explore the impact of modifications to the black 
hole production cross section, however, we consider cases both with and without this sup- 
pression. 
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FIG. 1: Cross sections a{vN -» BH) for n = 6 and (M D ,x min ) = (1 TeV,l) (solid), (1 TeV, 3) 
(long dash), (2 TeV, 1) (short dash), (2 TeV, 3) (dot-dash), and parton cross sections 7rr 2 (left) and 
irr'le~ I (right). The dotted curve is the SM cross section a(vN — > iX). The horizontal lines are 
the cross sections corresponding to interaction lengths equal to the vertical and horizontal depths 
of IceCube. (See text.) 



The neutrino-nucleon scattering cross section is then |LT | 



a(vN -»• BH) = 53 t 



where 



(a; min M D ) 2 /s 



M^/M D > 1 



dxai(xs) fi{x,Q) 



(6) 



(7) 



parameterizes the minimal black hole mass for which we expect a semi-classical description 
to be valid, s = 2m^E v , the sum is over all partons i, and fa are the parton distribution 
functions. The cross section is highly insensitive to choice of momentum transfer p"Tf ; 
e.g., the choices Q = M B h @ and Q = rj 1 lead to cross section differences of only 
10% to 20%. These cross sections are typically also insensitive to uncertainties at low x. 
For example, even extremely high energy neutrinos with E u ~ 10 10 GeV probe only x > 
(1 TeV) 2 /10 10 GeV 2 ps 10 -4 , within the range of validity of the CTEQ5 parton distribution 



functions we use 44 



Cross sections for black hole production by cosmic neutrinos are given in Fig. |T|. The SM 
cross section for uN — > £X is included for comparison (dotted curve). As a result of the sum 
over all partons and the lack of suppression from small perturbative couplings, the black 
hole cross section may exceed SM interaction rates by two or more orders of magnitude. 

Note that in our conventions, the cross section rises for increasing n and fixed Md- In 
conventions where the fundamental Planck scale is taken to be M* with 



1 

~G~d 



M. 



2+n 



(8) 



this behavior is reversed: the cross section decreases for increasing n and fixed M* |TT] . 
The dependence of the cross section on n for fixed Planck scale is convention-dependent 
and unphysical; we have adopted the convention to simplify comparison with existing 
collider bounds. 
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Assuming a constant density of 1.0 g/cm 3 for the Earth's surface near the detector, as is 
valid for IceCube, the neutrino's interaction length in Earth is 



1.7 x 10 7 km ( ^ J . (9) 



a 



The center of IceCube is at a depth of roughly 1.8 km. A neutrino reaching this point 
horizontally passes through 150 km of Earth. The cross sections corresponding to neutrino 
interaction lengths equal to these two lengths, that is, the horizontal and vertical depths of 
IceCube, are also given in Fig. |l|. We see that for the geometric cross section, Mb ~ 1 TeV, 
and neutrino energies E v ~ 10 9 GeV where the cosmogenic flux peaks (see Sec. [Hip , black 
hole production increases the probability of conversion in down-going neutrinos without 
increasing the cross section so much that vertical neutrinos are shadowed by the Earth. We 
therefore expect significantly enhanced rates in neutrino telescopes. (On the other hand, 
the upgoing event rates will be even more suppressed.) For the exponentially suppressed 
cross section, a similar enhancement is also possible for = 1 to 2 TeV. 

Once produced, these tiny black holes evaporate with rest lifetime of order TeV -1 ~ 
10~ 27 s. Even though highly boosted, they decay before accreting matter. They evapo- 
rate in a thermal distribution with temperature Tjj = (1 + n)/(47rr\.) |TD|, and average 
multiplicity M, H 



W - Mbh 2 * 
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Neglecting particle masses, the decay products are distributed according to the number of 
degrees of freedom 0: quarks (72), gluons (16), charged leptons (12), neutrinos (6), W and 
Z bosons (9), photons (2), Higgs bosons (1), and gravitons (2). We neglect the possibility 
of other low mass degrees of freedom, such as right-handed neutrinos and supersymmetric 
particles. About 75% of the black hole's energy is radiated in hadronic degrees of freedom, 
while the probability of any given decay particle being a muon (or a tau) is approximately 
3%. 

As through-going muons and taus will be a promising signal in neutrino telescopes, the 
typical multiplicity of black hole decays is of great importance. To quantify this, we define 
a weighted multiplicity 

A y =-E/ 1 dx(N)&i(x8)fi(x,Q) , (11) 

C j " , (ZminM D ) 2 /s 

where a is given in Eq. @, and (N) is as in Eq. (|ToD with Mg H = xs. The weighted 
multiplicity N is given in Fig. for various cases with n = 6 and = 1 TeV. For the 
geometric black hole cross section, these multiplicities may be substantially enhanced for 
lower n; for the exponentially suppressed cross section, the dependence on n is slight. We 
find that N ~ (9(10) is possible for ultra-high energy neutrinos. Note that while raising x m i n 
and including the exponential suppression factor both suppress the total cross section, they 
have opposite effects on iV: raising x m { n eliminates events with relatively low multiplicity, 
and so raises N, while the exponential suppression is largest for events with large M B h and 
large multiplicity, and so suppresses N. 
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FIG. 2: Weighted multiplicities N for n = 6, Mo = 1 TeV, and x x 
and parton cross sections 7rr 2 (left) and nr^e -1 (right). 
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III. ASTROPHYSICAL SOURCES OF ULTRA-HIGH ENERGY NEUTRINOS 



Observations of black hole formation and decay are most easily made with neutrino fluxes 
that are large near or above the EeV scale. Several such sources have been proposed in the 
literature. While a detailed discussion of these sources is beyond the scope of this paper, we 
briefly describe here a few of them, including those used in our study. 

First, neutrinos are almost certainly produced through pion decays in the scattering of 
protons off the cosmic microwave background, P7cmb — * niT+ |45|- This cosmogenic flux is 
subject to a number of quantitative uncertainties, including cosmological source evolution. 
As representative fluxes, we consider the recent results presented in Fig. 4 of Ref. These 
fluxes are shown in Fig. [| 

Second, gamma ray bursts have also been considered as a possible source of the highest 
energy cosmic rays. If this is the case, Fermi accelerated protons from shocks will generate 
extremely high energy neutrinos with energy spectrum d§ v /dE v oc E~ 2 p7j 48, |50 . 
This neutrino flux, as well as those from other compact sources, such as active galactic 



nuclei, is limited by the Waxman-Bahcall (WB) bound |51]. This constraint is valid for all 
astrophysical neutrino sources that are optically thin to pj and pp interactions. We consider 
a conservative estimate of this bound, E% d^ v jdE v = 1 x 10~ 8 cm -2 s _1 sr _1 [51|, where here 
v = v^v^Vp. This flux is also shown in Fig. |3|. It is approximately equal to the cosmogenic 



flux for E v ~ 10 10 GeV, but is much larger for lower (and higher) energies. 

Third, if the highest energy cosmic rays observed are generated by the annihilation of 
superheavy dark matter particles or from the decay of topological defects, neutrinos will also 
be produced. 



Such fluxes have been described in Refs. [|2|, |3|, pq| . We will not discuss 
these sources further, but note that they may predict large neutrino fluxes at extremely high 
energies, enhancing the results given below. 

In propagating to the Earth, the neutrino fluxes of Fig. 



will mix 56 . Given the solutions 



preferred by neutrino oscillation experiments and the enormous distances traveled, we take 
the neutrinos that reach the Earth to be in the ratio v e : : v T = 1 : 1 : 1 and similarly for 
anti-neutrinos. 
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FIG. 3: Representative fluxes: cosmogenic v (dashed) and v (dotted) |^6| and the Waxman-Bahcall 
flux *„ e = ^ = ^ (solid) p. 

IV. EVENT SIMULATION 

With regard to signals at neutrino telescopes, black hole evaporation products may be 
divided into three categories: showers (hadronic and electromagnetic), muons, and taus. 
Black holes also decay to neutrinos, but we neglect this flux in this work. 

The ranges of typical hadronic and electromagnetic showers are much less than the linear 
dimension of large-scale neutrino telescopes, and so, to first approximation, only contained 
showers from black holes produced inside the detector may be detected. The backgrounds for 
showers from black hole evaporation consist of hadronic showers from neutral and charged 
current neutrino events and electromagnetic showers from charged current electron neutrino 
events. At IceCube, hadronic shower energies should be measurable with an accuracy of 
about 30%. 

In contrast to showers, at the typical energies of black hole events, muons travel several 
kilometers before losing a decade of energy. The dominant signal is therefore through-going 
muons. IceCube can measure the energy and direction of any observed muon. The angular 
resolution is about 2° — 3° while the energy resolution is approximately a factor of three. 
Signal and background muons may therefore be differentiated with an energy cutoff. As we 
will see, for down-going muons with energy above 500 TeV to 1 PeV, the black hole signal 
may be well above the SM background from atmospheric neutrinos. Note that for black 
holes produced sufficiently near the detector, muon events may be obscured by showers. 
This occurrence is rare, however, and we ignore this possibility below. 

At high energies, when tau decay is sufficiently time dilated, taus have ranges as large 
or larger than muons, and so the dominant tau signal is from through-going taus. These 
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events have a characteristic signature consisting of a "clean" track, i.e., a track without 
much energy lost through low energy cascades. We will consider all events in which a tau 
track passes through the detector. Taus can be differentiated from muons by the cleanliness 
of their tracks or through their decays in the detector leading to "lollipop" events |57 



described later in this paper. We assume that taus and muons are distinguishable at all 
energies, although it will be difficult to distinguish a slow muon of energy < 200 GeV from a 
very high energy tau, as they can both generate clean, through-going tracks in the detector. 

Very massive black holes with large multiplicity decays may evaporate to several muons 
(or taus). However, these travel in coincidence through the detector. Their angular spread 
is only A6 ~ Th(N)/E u , much less than a degree, and so cannot be resolved with large 
Cerenkov detectors such as IceCube. Spectacular multi-lepton signatures are therefore not 
possible. Note, however, that more massive black holes are easier to detect, as they produce 
leptons with greater energy, which travel further before dropping below the cutoff energy. 

To evaluate black hole detection prospects, it is essential to determine, in a unified frame- 
work, both the SM background and the black hole event rate. We now describe both of these 
calculations. 



A. Standard Model Events 



1. Showers 

In the SM, a general expression for the total number of shower events in an underground 
detector is 



/f 
dcos6 z J dE Vi -^(E Ui ) Pi 



x / . . dy ) - ^{E Vi ) P mt , (12) 

where 9 Z is the zenith angle (8 Z = is vertically downward), and the sums are over neutrino 
(and anti-neutrino) flavors % = e,fi,T and interactions j = CC (charged current) and NC 
(neutral current). A is the detector's cross sectional area with respect to the v flux, T is its 
observation time, and dQ u JdE v . is the differential neutrino flux that reaches the Earth. For 
i = r, Eq. (|l"2"D is modified to include the effects of regeneration, as discussed below. 
P surv is the probability that a neutrino survives to reach the detector. It is given by 

P surv = exp[-X(e z )a^(E Ui )N A ] , (13) 

where N A — 6.022 x 10 23 g" 1 , and the total neutrino interaction cross section is 

= ^ cc + ^ NC ■ (14) 

Note that this is conservative, as it neglects the possibility of a neutrino interacting through 
a NC interaction and continuing on to create a contained shower. X(8 Z ) is the column 
density of material the neutrino must traverse to reach the detector with zenith angle 9 Z . It 
depends on the depth of the detector and is given by 

X{9 Z ) = f p{r{6 z ,l))dl , (15) 
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the path length along direction 6 Z weighted by the Earth's density p at distance r from the 
Earth's center. For the Earth's density profile, we adopt the piecewise continuous density 
function p(r) of the Preliminary Earth Model pB . 

P int is the probability that the neutrino interacts in the detector. It is given by 



P 



int 



exp 



. L J SM (E U .)_ 



(16) 



is the mean free 
L < L 3 m (E Vi ), and so P int « 



where, for showers, L is the linear dimension of the detector, and Lg M (P, 
path for neutrino interaction of type j. For realistic detectors 
L/ L J SM (E Ui ). To an excellent approximation, then, the shower event rate scales linearly with 
detector volume V = AL, and we present results in units of events/volume/time. 

Finally, the inelasticity parameter y is the fraction of the neutrino energy carried away 
in hadrons. The limits of integration are determined by the interaction type and neutrino 
flavor. For NC v e interactions and all and v T interactions, y max = 1 and y m \ n = E^/E„, 
where E^ 1 is the threshold energy for shower detection. For CC v e interactions, the outgoing 
electron also showers, and so y max = 1 and y m i n = 0. 



2. Muons 



Energetic through-going muons are produced only by CC interactions. For a muon 
to be detected, it must reach the detector with energy above some threshold P* hr . The 
expression of Eq. fll2f) then also describes the number of muon events with P surv as before, 
but with 

R 

Pnt = 1 - exp - TCr ,^ N , (17) 

. — y)E v and final energy P 



where is the range of a muon with initial energy P M = 
We assume muons lose energy continuously according to 



thr 
A 4 ' 



-a 



where a = 2.0 MeV cm 2 /g and (3 



dX 

4.2 x 1(T 6 cm 2 /g 



The muon range is then 



a + pE^ 
a + (3Ef* j 



(18) 



(19) 



In this case, y max = 1 - P* hr / E v and y min = 0. 

The event rate for muons is significantly enhanced by the possibility of muons propagating 
from several km into the detector. Note, however, that this enhancement is ^-dependent: 
for nearly vertical down-going paths, the path length of the muon is limited by the amount 
of matter above the detector, not by the muon's range. This is taken into account explicitly 
in the simulations, and its effect will be evident in the results presented in Sec. [VI]. 

At extremely high neutrino energies, the approximate form P int w R p / Lg^(E v ), often 
presented in the literature, is less accurate than the one used here, and the difference may 
be significant in scenarios in which the neutrino cross section is enhanced with respect to 
the SM value, such as the ones we explore in this paper. Note also that in this case, our 
expression for P surv is conservative, as it demands that the neutrino survive all the way to 
the detector, neglecting the possibility that it may convert a significant distance from the 
detector and still produce a signal. 
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3. Taus 



Taus are produced only by CC v T interactions. This process differs significantly from the 
muon case, as tau neutrinos are regenerated by tau decay through v T — > r — > v T As 



a result, for tau neutrinos, CC and NC interactions do not deplete the v T flux, but serve 
only to soften the neutrino energy. We include this important effect by first performing 
a dedicated simulation that determines E Ut (E Ut ,9 z ), the average energy a v T has when it 
reaches the detector, as a function of its initial energy E Vt and zenith angle 9 Z . The tau 
event rate is then given by 

d<&„, , „ s fy*™* , 1 dcr£9 ,— 



N T = 2tcAT / dcos6 z / ^ T ^BJ / dy CC/ - , f U (E„ T 

R T {{l-y)E UT ) 



exp 



e((l-y)E UT -E^) , (20) 



where R T ((l—y)E UT ) is the range of the produced tau, but evaluated at the energy of the tau 
neutrino after regeneration. R T is given by Eq. ([19]) but now with f3 = 3.6 x 10~ 7 cm _2 /g |p9" 



The last factor takes into account the requirement that the tau track be long enough to be 
identified in the detector. We require Ef 1 * ~ 2.5 x 10 6 GeV so that the tau decay length 
is above 125 m, the string separation length in IceCube. It is not clear, at this time, 
whether through-going tau events will be separable from less energetic muon events. Those 
tau events that include one (lollipop events) or two (double bang events) showers in the 
detector volume will be identifiable, however. The rate of down-going lollipop events in 
a km 3 neutrino telescope is expected to be of the order of the rate of down-going shower 
events, probably slightly smaller. Double bang events will be mostly observed for neutrino 
energies in a limited range between roughly 10 and 100 PeV |57|]. Our results assume that all 



tau events can be distinguished from muon events, although this may be difficult to realize. 

As with muons, at very high energies taus can travel several kilometers before decaying 
or suffering significant energy loss. The enhancement to tau event rates from this effect is 
# z -dependent as discussed above for muons. 

B. Black Hole Events 

Black hole event rates may be determined with only minor modifications. For showers, 
the corresponding expression is 



/f dsP 
dcos6 z dE n -^{E Vi )P mrv 

x f dy l — ^(E Ui ) P mt P BH (^) • (21) 

Jo a B H{E Ui ) dy 

The survival and interaction probabilities are now 

P surv = exp[-X(^)a tot (P, i )iV A ] (22) 

L 



-Pint = 1 - exp 



(23) 
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where 

a«* = a™+o BK , (24) 

with ct B h the cross section for black hole production, and L B n(E u .) is the neutrino mean 
free path for black hole production. As in the SM, L is the linear dimension of the detector 
in the case of shower events. We assume da^n/dy oc S(y — 0.75), i.e., that 75% of the black 
hole energy is carried away by showers, and impose that the generated shower has energy 
above E^ r to be detected. Pbk(E u ) = 1 for shower events. 

In the case of muon and tau events, L is the corresponding range for the average muon 
or tau energy. We assume E^ T = E u /N, where iV is the weighted multiplicity discussed 
in Eq. (JTTJ) . The true E^ T distribution is essentially flat with endpoints and 2E U /N. 
Our simplification is valid except for lepton energies near threshold. In addition, spreading 
the lepton energy has two compensating effects, as some leptons below threshold become 
detectable, and some above threshold drop below threshold. We have checked that the error 
made is insignificant at the ~ 10% level. 

To account for the branching fraction to leptons, we include 



Pbr(E Ui ) = 1 - exp 



n(e u : 

30 



(25) 



the probability of obtaining at least one muon (or tau) in the decay of a black hole when 
the expected muon (or tau) multiplicity is A/30. 



V. RATES 



The integrated event rates for showers, muons, and taus are shown in Tables |, ||, and |Tl 



respectively. We consider IceCube, with a representative depth of 1.8 km. We give results 
for various n, Md, Xmin, and with and without exponential suppression in the parton cross 
section. Only down-going rates are presented. Up-going rates are, as expected, extremely 
suppressed in the presence of black hole production, as will be seen in Sec. |Vj. 

For showers, the geometric cross section a = itt^, and Md = 1 TeV, we find a few events 
per year in each channel for cosmogenic fluxes, and as many as tens of events per year for the 
WB flux. These event rates are far above SM background. For the WB flux, bounds from 
AGASA and Fly's Eye will imply limits above Md = 1 TeV, but even for Md = 2 TeV, we 
find reasonable rates. In this case, however, the SM background will be important. For the 
exponentially suppressed parton cross section, event rates are suppressed, but not drastically 
so. In fact, for x min = 1, ~ Oil) event per year is possible given the cosmogenic flux. For 
a^min = 3, the exponential suppression is large, and event rates are highly suppressed. 

Generally, event rates are on more solid footing for larger x m i n , where both the semi- 
classical production cross section and the assumption of thermal decay are more reliable. 
For showers, however, our event rate calculation relies essentially only on the requirement 
that black holes decay to visible showers and is insensitive to the details of evaporation. At 
the same time, while the production cross section for black holes (or their stringy Planck 



mass progenitors ||42|| ) is subject to significant quantum corrections for Mbh ~ Md> there 
is no reason to expect it to vanish or be greatly suppressed. For showers, then, we find the 
requirement x m \ n = 1 reasonable. 

Qualitatively similar conclusions apply for muons and taus. While these event rates are 
suppressed relative to shower rates by the branching ratio for black hole decay to leptons, 
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TABLE I: Event rate for down-going showers (in 2ir sr) with E^ 1 = 500 TeV in IceCube. We 
consider the Waxman-Bahcall [^] and cosmogenic [^6| fluxes, Md = 1 and 2 TeV, and various 
cases (n,x m [ n ,a), where n is the number of extra dimensions, x m [ n = Mg^/iVfo, and a is the 
parton level cross section for black hole production. 



Showers (km 3 yr 1 


WB Flux 


Cosmo^ 


;enic Flux 


M D = 1 TeV 


M D = 2 TeV 


M D = 1 TeV 


M D = 2 TeV 


Standard Model 


4.8 


4.8 


0.1 


0.1 


BH (6,l,vrr 2 ) 


44.6 


3.1 


5.2 


0.9 


BH (6,3,vrr 2 ) 


6.5 


0.5 


2.3 


0.3 


BH (6,l,7rr 2 e- 7 ) 


18.5 


1.2 


2.1 


0.3 


BH (6,3,7rr 2 e~ 7 ) 


0.4 


2.8 x 10~ 2 


0.1 


1.5 x 10~ 2 


BH (3,1, tit?) 


16.3 


1.1 


2.7 


0.4 


BH (3,3,^r 2 ) 


3.0 


0.2 


1.2 


0.1 


BH (3,l,7rr 2 e- 7 ) 


3.8 


0.2 


0.5 


6.1 x 10~ 2 


BH (3,3,7rr 2 e- 7 ) 


2.9 x 10" 2 


1.9 x 10~ 3 


8.6 x 10~ 3 


9.3 x 10" 4 



TABLE II: As in Table |, but for the flux of down-going muons and Ef r = 500 TeV. 



Muons (km 2 yr 1 ) 


WB Flux 


Cosmos 


;enic Flux 


M D = 1 TeV 


M D = 2 TeV 


M D = 1 TeV 


M D = 2 TeV 


Standard Model 


6.0 


6.0 


0.2 


0.2 


BH (6,l,vrr 2 ) 


27.7 


2.6 


4.9 


1.1 


BH (6,3,vrr 2 ) 


12.0 


1.1 


4.9 


0.7 


BH (6,l,7rr 2 e- 7 ) 


8.8 


0.7 


1.3 


0.2 


BH (6,3,7rr 2 e- 7 ) 


0.6 


4.6 x 10~ 2 


0.2 


2.5 x 10~ 2 


BH (3,1, tit?) 


14.0 


1.2 


4.2 


0.6 


BH (3,3,vrr 2 ) 


7.1 


0.6 


3.5 


0.4 


BH (3,l,7rr 2 e- J ) 


1.8 


0.1 


0.3 


4.1 x 10~ 2 


BH (3, 3, tit 2 e" 1 ) 


4.2 x 10~ 2 


3.1 x 10~ 3 


1.5 x 10" 2 


1.6 x 10~ 3 



TABLE III: As in Table 



EL but for the flux of down-going taus and £* hr = 2.5 x 10 6 GeV. 



Taus (km 2 yr , 



WB Flux 



Cosmogenic Flux 



M D = 1 TeV M D = 2 TeV M D = 1 TeV M D = 2 TeV 



Standard Model 



0.9 



0.9 



0.1 



0.1 



BH (6,l,7rr 2 ) 
BH (6,3,^r 2 ) 
BH (6,l,7rr 2 e- 7 ) 
BH (6,3,7rr 2 e- 7 ) 



15.1 
7.7 
4.5 
0.5 



2.2 
1.0 
0.6 

4.4 x 10~ 2 



5.0 
4.5 
1.3 
0.2 



1.2 
0.7 
0.2 
2.9 x 10~ 2 



BH (3,1, tit?) 

BH (3,3,vrr 2 ) 
BH (3,l,7rr s 2 e- 7 ) 
BH (3,3,7rr 2 e- 7 ) 



8.7 
4.8 
1.0 

2.9 x 10~ 2 



1.1 
0.6 
0.1 

2.8 x 10~ 3 



4.1 
3.1 
0.4 

1.6 x 10~ 2 



0.7 

0.4 
5.2 x 10" 2 
1.8 x 10" 3 
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they are enhanced by the possibility of muons and taus propagating from many km away 
into the detector. We find that these effects effectively balance each other. In the case of 
muons and taus, we have assumed thermal decay distributions, which may not be accurate 
for x min ~ 1. The event rates are hardly reduced for x min = 3, however; as noted in Sec. II, 
the more stringent x m [ n requirement preferentially eliminates events with low multiplicities, 
and so has little impact on lepton rates. 

The muon flux is significantly larger for the WB flux than for the cosmogenic flux. This 
is straightforward to understand — the WB flux is larger at low energies. The tau event 
rate is not as greatly enhanced, however. The tau decay length is 4.9 km(£' r /10 8 GeV). 
While the WB flux enhances fluxes at low energies, the resulting taus decay before they lose 
energy, and the tau range is therefore diminished. Larger low energy fluxes therefore do not 
enhance the tau rate as significantly. 

In addition to enhancing the SM event rates by more than an order of magnitude, Tables |- 
[IJ show that black hole production may also change the relative event rates of the various 
channels. The ability of neutrino telescopes, unique among cosmic ray experiments, to 
differentiate showers, muons, and possibly taus, allows one to measure these relative event 
rates. Given sufficient statistics, this may provide an important signal of physics beyond 
the SM. 



VI. ANGLE AND ENERGY DISTRIBUTIONS 

Neutrino telescopes may also measure angle and energy distributions, providing additional 
opportunities to distinguish black hole events from SM or other possible physics and for 
constraining black hole properties. 

The zenith angle distributions of SM and black hole events for showers, muons, and taus 
are given in Figs. [|, |5|, and |6], respectively. Black hole production makes the Earth even 
more opaque to ultra-high energy neutrinos, reducing the up-going rate, while increasing the 
probability of interaction, and thereby the event rate, of down-going neutrinos. In all cases, 
up-going events rates are below SM contributions. No sensitivity to black hole production 
is then expected when looking for up-going events. 

The characteristic features of the zenith angle distribution of showers compared to the 
distributions of muons and taus are also noteworthy. For showers, vertical fluxes are less 
attenuated than horizontal ones, and so the angular distribution is maximized for vertically 
down-going events. In contrast, as noted in Sec. [IV Aj, vertically down-going muon and tau 



events cannot benefit from muon and tau ranges beyond the depth of the detector (~ 2 
km). The optimal direction for muon and tau events is therefore closer to the horizon. In 
scenarios with enhanced cross sections such as the one we are exploring in this paper, the 
neutrino flux is attenuated even in the horizontal direction as can be seen in Figs. || and ||. 
As a consequence, the quasi-horizontal direction with small but positive cos# 2 , where the 
effects of attenuation are maximally offset by lepton range, leads to the largest rates. For 
muons, this optimal zenith angle is around cos #.2 ~ 0.2. For taus, with even longer ranges, 
it is cos6 z pa 0.4. These are general considerations that apply to the detection of ultra-high 
energy neutrinos in any scenario predicting enhanced cross sections with respect to the SM 
cross sections. 

Energy distributions for the various channels are given in Figs. [7|, and |9|. The en- 
ergy distributions are clearly sensitive to x min and the presence or absence of exponential 
suppression in the parton cross section. Different types of events therefore probe different 
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FIG. 4: Zenith angle distribution of shower events for 
(n,x min ,a) = (6, l,7rrf) (solid), (6,3,vrr;?) (long dash), (6,l,7rr^e 
(dot-dash). Also shown is the standard model prediction (dotted). 
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FIG. 5: As in Fig. |, but for muon events with £* hr = 500 TeV. 



and complementary aspects of black hole production and decay. By measuring the down- 
going energy distributions with reasonable statistics, the IceCube detector may be able to 
discriminate between the different possibilities. Note that to facilitate comparison with con- 
ventional presentations of SM rates at neutrino telescopes, we have plotted distributions in 
initial neutrino energy. Shower energies are related to these neutrino energies in a fairly di- 
rect way, and the shower energy distributions will have roughly the same shape. The energy 



14 



Si 



O 

o 

\ 
Z 
T3 




FIG. 6: As in Fig. |, but for tau events with E^ hr = 2.5 x 10 b GeV. 
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FIG. 7: Energy distribution of down-going shower events for Md = 1 TeV, EQ* = 500 TeV, and 
(n,x m - m ,a) = (6, l,vrr^) (solid), (6, 3,7rr;?) (long dash), (6, 1, ■nr 2 s e~ I ) (short dash), (6, 3, 7rr^e _/ ) 
(dot dash). Also shown is the standard model prediction (dotted). 



distributions of through-going leptons may be significantly distorted, however, as the lepton 
energy is a function of black hole multiplicity as well as the distance the lepton propagates 
before reaching the detector. Note also that, while shower energies should be well- measured, 
lepton energy measurements present significant challenges, especially for taus. Lollipop tau 
events will be a clear tau signature with well-measured energy. However, determining the 
energy of through-going tau events will be difficult in IceCube as they typically do not ra- 
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FIG. 8: As in Fig. fj], but for down-going muon events with = 500 TeV. 




diate and, therefore, appear similar to minimum ionizing muons with energy ~ 200 GeV. 
More detailed detector simulations will be necessary to effectively use the information from 
through-going taus described in this paper. 



VII. PROSPECTS AND CONCLUSIONS 

In the presence of TeV-scale gravity and extra dimensions, ultra-high energy cosmic 
neutrinos may produce microscopic black holes in the surface of the Earth. We have explored 
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the prospects for detecting such events in neutrino telescopes. We considered contained 
showers, through-going muons, and also through-going taus. 

The rates for a km 3 detector, such as IceCube, are of the order of a few per year in 
each channel for near current bounds. These rates are well above SM backgrounds, and 
provide a significant opportunity to observe black hole production. In several years of data 
taking at a kilometer-scale neutrino telescope, it should be possible to probe most of the cases 
described in this paper: values of Mp well above 1 TeV may be within reach, and observable 
event rates are even possible for exponentially suppressed parton cross sections. At the 
same time, these results imply that black hole observation at ANTARES and NESTOR 
appears unlikely, given their effective areas of ^ 0.1 km 2 . The RICE experiment ||61| , which 
aims to detect coherent Cerenkov radiation from electromagnetic and hadronic showers at 
frequencies of 100 MHz to 1 GHz, provides another interesting and complementary probe 
with effective volume for shower detection reaching ~ 1 km 3 above a few PeV. 

The expected rates for black hole events are, of course, highly sensitive to the choice of 
neutrino flux. As this work was being completed, a paper also discussing black holes at 
neutrino telescopes appeared p6 |. There, results for two fluxes are presented: a "limit from 
hidden sources" flux and a cosmogenic flux. The first flux is roughly 3 times higher than the 



limit placed by the AMANDA experiment on diffuse fluxes [B2|. The WB flux we use is much 



more conservative and is roughly 2 orders of magnitude smaller. A variety of cosmogenic 



fluxes 63, Em are also considered in Ref. 26 . Our ratios of black hole to SM rates are 



in agreement, but here, too, it appears we have chosen a more conservative representative 
flux, leading to a factor of 3 to 5 fewer shower and muon events from cosmogenic neutrinos. 
The possibility of tau events, exponentially suppressed cross sections, and angle and energy 
distributions were not addressed in Ref. J2B] . 

Relative to IceCube, the prospects for black hole detection appear to be slightly brighter 
at the Auger Observatory, expected to be completed by 2004, where tens of events per year 
may be discovered [TTL [TJ, OJ. IceCube is scheduled to reach its ultimate goal of 1 km 3 in 
2009, but it will collect data continuously as it grows in stages, beginning with the existing 
AMANDA II detector, roughly of area 0.1 km 2 , and reaching a volume of around 0.5 km 3 in 
2006. It is difficult to model this time-varying detection volume, but by 2007, an integrated 
exposure of roughly 1 km 3 yr will have been achieved. Given the rates of Tables fl-pl, we 
conclude that, for scenarios with Mjj near current bounds, IceCube may detect black hole 
events along with the Auger Observatory before the LHC begins operation. 

Of course, enhanced event rates alone are not necessarily a signal of new physics. In 
particular, the uncertainty of astrophysical neutrino fluxes makes it difficult to determine 
whether observed events result from the SM or new physics based on counting experiments 
alone. Comparisons of quasi-horizontal showers and Earth-skimming neutrino rates provide 



a powerful discriminant in ground arrays [11, 15]. However, the possibility of measuring 



relative event rates in three separate channels at IceCube provides a direct, unique, and 
complementary tool for distinguishing black hole events from the SM and other new physics 
possibilities. As noted in Sec. [VI], angular and energy distributions will also be helpful 
to separate the signal of a large astrophysical flux from the signal of new physics. The 
separation of these channels in ground arrays and fluorescence detectors is extremely difficult. 

Given sufficient statistics, then, the information provided by IceCube may provide in- 
formation on black hole branching ratios, and the various angle and energy distributions 
may also help distinguish various properties of black hole production and decay. Black hole 
production may be the best prospect for discovering extra dimensions in the near future, 
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and neutrino telescopes may be an excellent tool for studying them. 
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